ELSEVIER 


Journal of Power Sources 166 (2007) 400-410 


Available online at www.sciencedirect.com 


ScienceDirect 


JOURNAL OF 


www.elsevier.com /locate /jpowsour 


Nanoparticle gasifier fuel cell for sustainable energy future 


Dimitri Gidaspow*, Veeraya Jiradilok 


Illinois Institute of Technology, Chicago, IL, United States 


Received 6 December 2006; received in revised form 20 December 2006; accepted 21 December 2006 
Available online 17 January 2007 


Abstract 


A new concept for production of electricity from biomass or coal using molten carbonate fuel cells is proposed. It involves feeding fine coal 
particles or biomass, for sustainable energy future, with steam into the anode compartment of the fuel cell in which the waste heat from the fuel 
cell is used to produce synthesis gas which reacts electrochemically. This concept is illustrated using carbon nanoparticles as the fuel. 

A computational fluid dynamics (CFD) model for the gasifier-fuel cell has been developed. Concentration, temperature and current density 
profiles have been computed. The computations show that practical current densities can be achieved. However, for the new concept to work 
well the carbon monoxide produced by the gasification should be allowed to react electrochemically. The computed temperature distributions for 
adiabatic operation show an initial drop in temperature due to gasification, followed by a rise which will have to be balanced by staggering the 


cells or by other means using the CFD design method. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


There is an urgent worldwide need to develop sustainable 
energy sources and means of production of nonpolluting elec- 
tricity. Biomass is one such source of energy. It has been known 
for over a century than biomass particles, such as saw dust, can 
be used a source of fuel in a high temperature fuel cell to pro- 
duce electricity. Unfortunately a direct feed of particles into a 
fuel cell gives very low current densities. Fig. 2.7-2 in Liebhaf- 
sky and Cairns’s book [1] shows that fuel particles fed into a fuel 
cell do produce appreciable current. But for practical applica- 
tions the current density must be increased by almost two orders 
of magnitude. This should be possible by feeding biomass or 
coal nanoparticles with steam into a high temperature fuel cell 
(Gidaspow’s Patent application [2]). 


2. Biomass gasification integrated with fuel cell 
literature 


The use of carbon nanoparticles has recently been sug- 
gested [3-6] for the development of a direct carbon fuel cell, 
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C+Or & CO,, to take advantage of the fact that in an ideal 
situation all of the enthalpy of combustion of carbon can be 
transferred into electrical work in the reversible cell. In this 
innovative concept the carbon is dissolved in the molten car- 
bonate electrolyte which is recirculated in the fuel cell and is 
oxidized at the anode. Although the concept proposed here is 
different, the idea of circulating the molten carbonate may be 
used to remove impurities, such as ash and sulfur found in coal. 
The solids produced will be removed using steam with a cyclic 
operation to clean the electrodes. The sulfur and ash will then 
removed by recirculating and filtering the electrolyte. 

U.S. Patent 6,680,137 B2 [7] described a biomass gasification 
and fuel cell system. The fuel gas from the gasifier is directed 
to the anode of the fuel cell and at least a portion of the exhaust 
gas from the anode is directed to the combustor. The idea of 
gasification inside the fuel cell compartment is not described. 
This important new idea proposed here eliminates the expensive 
cryogenic air separation system in the Department of Energy 
FutureGen System for producing hydrogen from coal and using 
hydrogen in a fuel cell. 

Donolo et al. [8] and Tomasi et al. [9] described the use of 
biomass in molten carbonate fuel cells. The gases are produced 
in separate gasifiers and reformers. Panopoulos et al. [10,11] 
described the use of a heat pipe to transfer heat from the fuel 
cell to the gasifier. Morita et al. [12] described electrochemical 
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Nomenclature 

Ca drag coefficient 

Ck heat capacity of phase “k” (Jkg~! K7!) 

Cp heat capacity of gas (J kg7™! K7!) 

C pi initial concentration of fixed carbon in the parti- 
cles 

dp particle diameter (m) 

Di bulk gas diffusivity of gas species i (m? s7!) 

De.i effective diffusivity of gas species i in the core 
(m? s7!) 

DmM,i effective diffusivity of gas species i in the shell 
(ui? s7!) 

E(P;) reversible emf of fuel cell (V) 

Eo Standard emf of fuel cell (V) 

F Faraday constant (C kmol” !) 

g gravitational acceleration (m s7?) 

go radial distribution function 

Gs solid modulus (Pa) 

h heat transfer coefficient (kg s73 K7! m7!) 

hyk heat transfer coefficient (kg s73 K7!) 

H enthalpy (Jkg~!) 

I current density (A m~?) 

Kg gas phase heat transfer (kg ms~* K7!) 

Kp,i mass transfer coefficient of gas species i through 
gas film surrounding a particle 

K; rate constant (kmol (kmol C)~! atm~! s7!) 

Ks particulate phase heat transfer (kg m s7? K7!) 

m* mass production or consumption (kg m~? s7!) 

Mi molecular weight of species i (kg k mol~!) 

n number of electrons 

Nu Nusselt number 

P continuous phase pressure (atm) 

Pr Prandtl number 

Pi partial pressure (atm) 

P; particulate phase pressure (atm) 

P“ equilibrium pressure (atm) 

R gas constant (8.314 kJ kmol~! K7!) 

Re Reynolds number 

Reff effective cell resistance of a molten carbonate fuel 
cell (ohm m?) 

Ri rate of reaction of reaction i (kmol m~? s7!) 

t time (s) 

T temperature (°C) 

V velocity (m s7!) 

Vm terminal cell potential (V) 

We weight fraction of carbon 

Y; weight fraction of species i 

Xi mole fraction of species i 


Greek letters 


anode channel thickness 

interphase transfer coefficient 
stoichiometric coefficient of reaction i 
volume fraction of phase “k” 

effectiveness factor for the reaction in core 


Oc void fraction of the particles in core 
Os void fraction of the particles in shell 
u viscosity 

p density 

PUC fraction radius of the unreacted core 
Qi Thiele modulus for reaction i 
Subscripts 

s solid phase 

g gas phase 


production of current in a molten carbonate fuel cell with oxida- 
tion of hydrogen and carbon monoxide produced from biomass. 
At IIT Chuck [13] using a molten carbonate cell manufactured 
at IGT (Institute of Gas Technology), now GTI (Gas Technology 
Institute), systematically measured the current density of oxida- 
tion of carbon monoxide as a function of polarization from 10 to 
70% inlet carbon monoxide concentration and temperatures of 
750-875°C. All the data were correlated by the relation current 
density equals Nernst potential minus voltage divided by cell 
resistance. The same model is used here. 


3. Inherent efficiency of molten carbonate fuel cell 
systems for electricity production 


Production of electricity from coal, biomass or organic waste 
according to the second law thermodynamics is best done using a 
fuel cell. In a fuel cell the differences of chemical potential of the 
fuel and air are forced to produce electric work [1]. In combus- 
tion the fuel and air are mixed producing heat, a random form of 
energy. Then this randomness is used to produce work using vari- 
ous cycles. Hence, inherently such systems are less efficient than 
those utilizing fuel cells. This is the advantage of gasification of 
fuels over combustion, as in the FutureGen project [14]. Unfor- 
tunately the present version of FutureGen involves gasification 
with oxygen to enable sequestration of carbon dioxide. 

The new concept proposed here eliminates the oxygen plant 
and uses the waste heat from the fuel cell to supply the energy for 
production of synthesis gas in the same unit eliminating entropic 
losses due to heating and cooling. The molten carbonate fuel 
cells require carbon dioxide in the cathode for the transfer of 
carbonate ions [1]. Hence, they are the natural fuel cells to be 
used with fuels containing carbon, such as natural gas, coal or 
biomass. When the carbon dioxide is mixed with air in the cath- 
ode compartment at high concentrations, it provides additional 
electricity due to the differences in concentrations in the anode 
and the cathode. The carbon dioxide can be transferred from the 
anode to the cathode using a mass-exchanger in which the carbon 
dioxide reacts with potassium carbonate contained in a porous 
Teflon sheet [15]. The carbon dioxide is produced by the ther- 
mal decomposition of the product, bicarbonate using waste heat. 
Furthermore the carbon monoxide that will be present during the 
reaction with any carbonaceous fuel is itself electrochemically 
oxidized, and does not act as a poison [13]. 
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Fig. 1. The large random kinetic energy, Granular temperature, of nanoparticles 
[18] compared to micron size particles [39-41]. 


4. New concept 


The new concept is based on our recent discovery of the 
unique flow properties of nanoparticles [16—18]. Nanoparticles 
fluidize without formation of bubbles due to their Brownian 
type interaction with air. Fig. 1 shows that their random kinetic 
energy, called granular temperature, is very high at rather low 
gas velocities. This random motion, of the order of 1 m? s7? pro- 
duces a mass flux equivalent to a current of 7 A cm7?. Hence, 
mass transport of particles will not limit the performance of the 
fuel cell. The rate of gasification of the nanoparticles will also 
be very high due their small size. 

Hence, the proposal is to feed carbon or biomass particles 
with steam into a molten carbonate fuel cell operating at about 
600-800 °C, as shown in Fig. 2. In the fuel compartment of the 
fuel cell water will be split using biomass or carbon nanoparti- 


Biomass 

or 
Carbon 
nanoparticles 


Steam 


ANODE 
CURRENT COLLECTOR 


cles, with the additional energy supplied by the heat generated 
by the fuel cell. The hydrogen generated will react electro- 
chemically to produce current. Unreacted fuel will be burned to 
maintain the high temperature and to provide the carbon dioxide 
to the cathode of the molten carbonate fuel cell. 

The fuel cell will be initially similar to the molten carbonate 
fuel cell being commercialized by the Fuel Cell Energy Corpora- 
tion, founded by the late Baker [19]. In these fuel cells, natural 
gas, essentially methane is internally reformed with steam to 
produce hydrogen which then reacts at current densities of over 
100 mA cm~? and an efficiency of about 50%. The concept has 
been scaled up to 1.8 MW [20] and is near commercialization. 
The novel idea proposed here is to use a storable fuel, carbon or 
biomass nanoparticles. Such a fuel cell battery can be used to 
power an automobile or a tractor without the need of an addi- 
tional fuel. Nanoparticles will be stored in the fuel tank and 
delivered to the fuel cell with steam or carbon dioxide. 

Commercial carbon nanoparticles are available. Biomass 
fuel, such as switch grass, has a high reactivity for reaction with 
steam [21]. 


5. Computational fluid dynamics (CFD) model 
5.1. Hydrodynamics model 


The CFD model is a predictive hydrodynamics multi-phase 
model developed to model fluidization [22,23] and nanoparticle 
flow [17]. In view of the measured nanoparticle viscosity and 
stresses we used the viscous model in Jiradilok et al. [17] to 
predict the hydrodynamics of nanoparticles in a gasifier fuel 
cell. The equations are the conservation of mass, momentum and 
energy equations for each phase given in Gidaspow’s book as 
model B. All equations are written in rectangular coordinates, as 
summarized in Table 1 . These are 4(N + 1) nonlinear-coupled 
partial differential equations for 4(N + 1) dependent variables, 
for | gas phase and N solids phases. The variables to be computed 


Unreacted fuel to 
steam generator and 
combustor 


MOLTEN SALT 


ELECTROLYTE, 600 — 800°C 


Burned anode gases 


Preheated air 


Current 


CATHODE: 
z +CO, +2e & CO;- 


Heat Exchanger 


Fig. 2. Nanoparticle gasifier fuel cell. Anode reactions [34,33], gasification (for carbon feed), carbon steam reaction: C+H20— CO+Hp) (endothermic 
AH=32 kcal gmol~!); water shift reaction: CO + H20 <> CO2 + H? (exothermic AH = —9.8 kcal gmol—!); electrochemical reaction: H? + co,” > H20 + CO? 4 


2e (exothermic —AH=(E — VJI — (T AS/nF)I. 
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Table 1 
Mathematical model 


Continuity equations 


0 : 
ge +V: (pi£ivi) =m; (i = gas or solid) 


Momentum equations 

O(Pg&g Ug) 
ðt 

VPI +V -Tg — PB(Vg — Vs) + Peg 

ƏlPsEsVs) 
or 


+ V  (PgEgVgvg) = 


+ V + (Ps&sUsUs) = 
N 


VRI+V-% + pes- vs) += | o Y erpe) g 


E 
g 
k=g,s 
Energy equations 


N 
È (Ep Pp Ha) + V+ (€gPe Hve) = (42 + vg- VP) + XOA- 
k=1 


H>,CO 


Tg) +V. (Kg£gVTa) + eS Rig AHig T eX A Hi fuel cell 


i i 


a 
g EPs Hs) F V. (Es ps Hg vs) = h(Tg = Ts) + V. (Kse,VT;) 


Constitutive equations 
(1) Definitions: eg + £s = 1 


(2) Gas pressure: Pg = PeRTg 


(3) Stress tensor (i= gas or solid) 


= = 2 =A 1 
Ti = 2D; + (u — Zu) t(D) with Di = 5 [Vv + (Vv;)"] 


(4) Empirical particulate phase viscosity and stress model 


V P; = Glég)Ves with G(eg) = 107!4926£8+18-667 dyn cm~? 


ps = 0.01782! go (Poise) 


(5) Fluid-particulate interphase drag coefficients 
for e <0.8 (based on the Ergun equation) 


2 

Ey Es 
p= 150548 4 1,7525 

egds Egdp 


for € > 0.8 (based on the empirical correlation) 


3 
B= 7e 


[vg Us| 


Pgés|Ug — Us| 972-65 
g 
dp 


24 0.697 
where Ca = Re + 0.15Re5 ] for Rep < 1000 and Cg = 0.44 
e 


for Rep > 1000 


(6) Enthalpy 
Hz = Cp, (Tz — Ty) 


H, = Cp, (Ts — Tp) 
where Ty T? : the standard temperature, 25°C, 298 K 


(7) Heat transfer coefficients 
(7.1) Gas-particle heat transfer coefficient, hyk 

for s < 0.8 

Nux = (2 + 1.1Re®® Pr!/3) Sp (Re < 200) 
4Re 0.83 

= 0.123 (=) 50-17 (200 < Re < 2000) 
= 0.61 ROST g (Re > 2000) 

for £ > 0.8 
Nux = (2 + 0.16 Re’), (Re < 200) 
= 8.2Re”6 Sp (200 < Re < 1000) 


= 1,06Re®-457 5, (Re > 1000) 
Eg Pgl¥g—Vg dx Pr= Cp.gHg 
Hg i K, 


Nu = "tt, h= Shu [36] 


where Re = 


6 
; Sk = Ek T> 


Table 1 (Continued ) 


(7.2) Gas heat transfer coefficient 


au 1.786 
0 _ 3-1 
Ky = 8.64 x 10 (z) (gcms °K) 


Ke = (1 — Je)Ko 


(7.3) Particle heat transfer coefficient, Ky [23 
Kx re 
== Va) + ve [oR +O - 0 =| 
* 


g 


with “2 = 
À 


2 or et B! +) : 


1-B/R\ G- B/R} °R 1-B/⁄R 2 


K 1 10/9 

=e 
>, a= 125( :) , G=7.26 x 107° 
K: Eg 


g 


are the pressure P, the solids volume fractions ex (k= 1, .. ., N), 
the gas velocity components Ug and Vg and the solids velocity 
components U% and Vz (k=1,...,N) in the —x and —y direction, 
respectively, and the temperature Tg and T; (k=1,...,.N) of gas 
and solids phases, respectively. The numerical scheme used in 
the IIT code is the Implicit Continuous Eulerian (ICE) approach. 
The model uses donor cell differencing. The conservation of 
momentum and energy equations are in mixed implicit form. The 
continuity equations excluding mass generation are in implicit 
form. In this study there is one solids phase. 

The solid viscosity and solids modulus are input data into 
the viscous model. The solid viscosity increases with increasing 
solid concentration, which was estimated as a function of radial 
distribution function [17]: 


us = 0.0017e!/3 ¢9 (Pas) (1) 


The above expression for the solids viscosity is based on the 
experimental data. The radial distribution function at contact, 
go, is calculated as follows: 


‘ 1/3] 7! 
go = |1-|— (2) 
Es max 


where the maximum solids packing, £s,max Was estimated from 
the solid compression in the IIT fluidized bed and simulation 
data which was 0.08. 

The correlation for particulate viscosity given by Eq. (1) is a 
semi-empirical equation based on kinetic theory. The viscosity 
increases to the one-third power due to isentropic compression 
[24]. The radial distribution function is based on Bagnold’s equa- 
tion. It is similar to the viscosity of a concentrated suspension 
of Frankel and Acrivos [25]. 

The solid stress modulus of nano-size particle for 10 nm Tul- 
lanox was estimated by Jung and Gidaspow [26]. They measured 
the solid volume fraction as a function of bed height using a y-ray 
densitometer in a settling experiment giving the solid pressure 
as a function of the solids volume fraction. Hence, the deriva- 
tive of solid pressure equals the solid stress modulus. The solid 
modulus is expressed as a function of the gas volume fraction. 
An exponential form is used. 


G(ég) = 107 14.926, +18.667 (dyn cm~?) (3) 
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The gas and solids phase energy equations are written in their 
enthalpy form. The viscous dissipation of two phases and the 
slip work between phases were assumed to be not important 
and were not included in the energy equation. There are two 
important heat sources, endothermic heat from gasification reac- 
tions and exothermic heat from fuel cell. Therefore, in the fuel 
compartment of the fuel cell water will be split using carbon 
nanoparticles, with the additional energy supplied by the heat 
generated by the fuel cell. The hydrogen and carbon monoxide 
generated will react electrochemically to produce current. 


5.2. Gasification 


5.2.1. Description of gasifier 

The gasification reactions are similar to those for a moving- 
bed coal gasification reactor of Yoon et al. [27]. The dominant 
processes of the gasification are heterogeneous reactions 
between carbon and components of the gas phase. 

The gasification reactions consist of three reactions as fol- 
lows: 


Reaction! : C + CO2—~ 2CO 
Reaction2: C + 2H —> CH4 
Reaction3: C + H20 —> CO + H2 


In addition to the three heterogeneous reactions, the water 
gas shift reaction occurs in the gas phase catalyzed by carbon 
particles. 


Water shift reaction: CO + H20 & CO2+H2 


Skinner and Smoot [28] have reviewed the literature of the reac- 
tions of char and carbon. 


5.2.2. Model 

For the gasification the conservation of species equations are 
added into the code to describe gasification reactions. The left- 
hand side terms of the continuity equations for solids and gas 
phase do not equal to 0 due to mass changing with reactions. 
They are summarized as follows: 


Continuity equation for solid phase: 


0 
gee» + V - (psésVs) = Me(—Rı — Ro — R3) (4) 


Continuity equation for gas phase: 


“(ost + V : (PgEg Vg) 
= Rı(—Mco, + 2Mco) + R2(—2My, + Mcuy) 
+R3(Mco + Mn, — Mu,0) 
+R4(Mco, — Mco + Mp, + Mp,0) (5) 


e The conservation of species balance in solid phase: 


component balance for C: 


ð 
geses We) + V - (Ps€s Wc Vs) = Mce(—Rı — R2 — R3) 


(6) 
e The conservation of species balance in gas phase: 
component balance for CO? : 
ð 
g Vefe Yco) +V. (PgEgYco, Va) = Mco, (—Rı + R4) 
(7) 


component balance for CO: 


0 
ay Paee¥co) + V - (Pg£€gYco Vg) = Mco(2Ri + R3 — R4) 
(8) 


component balance for H? : 


a 
g (Patti) + V : (Og£gYn, Vg) = Mm, (—2Ro + R3 + Ra) 
(9) 


component balance for CH4 : 

ð 

g PefeYcm) + V - (PgE£gYcH, Ve) = Mcu, R2 (10) 
component balance for H20 : 


ð 
g Pe£e Ymo) + V > (Pg£€g Yn o Vg) = Mpm,o(— R3 + R4) 
(11) 


where Y; is the weight fraction of species i in gas phase and 
Wc is the weight fraction of carbon. 


5.2.3. Kinetic study 
5.2.3.1. Heterogeneous reaction model. The shrinking core 
model used to calculate the rate of the heterogeneous reaction 
is given by: 
A es(Pi — PP) 

o (dp/6Kp,i) + 31 — puc)RT/12puc Dm,i) + (1/nipicKriCe) 


Similar to that given by Yoon et al. [27] there are three main 
resistances in this model, mass transfer, diffusion and reaction 
effects. In the system of flow of nanoparticles, the first two 
effects, mass transfer and diffusion are small due to the small 
particle diameter. The reaction part in the shrinking core model 
plays the important role. 

The effectiveness factor for the reaction in the core, 7j;, is 
defined by: 


4 l 1 E 
S (am a) ) 


The Thiele modulus for reaction i is defined by: 


d K,iC 
ĝi = pPUC ric (14) 
6 Yi De,i/ RT 


(12) 
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The effective diffusivities in the outer shell and also in the core 
of the particles are estimated by Walker et al. [29] as follows: 


Shell: Dm, = Di6§ (15) 
Core: De; = D;0% (16) 


The Gasification reactions consist of three heterogeneous reac- 
tions as follows: 


Reaction 1: C+CO, —> 2CO 


(AHR 1 |298k = 40.273 kcal gmol') (17) 
45, 000 
K;.co, = 930 exp (- RT; ) 
(gmol (gmol C)~! atm™! s7!) (18) 
P2 40, 300 
Keg = -®© = 1.222 x 10” exp (- ) (19) 
Po, RTs 
P2 
Pe co (20) 


1.222 x 10° exp (—40, 300/ RTs) 


The rate constant, K;,co,, is for Pittsburgh 8 coal [30]. Skinner 
and Smoot [28], Table 2 gives the Arrhenius constants for reac- 
tion 1 for charcoal, coke etc. The equilibrium constant, Keg, is 
the one given in Yoon et al. [27]. 


Reaction 2: C+ 2H —> CH4 


(AAR 3l298K = —21.849 kcal gmol”') (21) 
E 1650 

K.m, = 8.36 x 10 äi exp (5) 

(gmol (gmol C)~! atm™! s7!) (22) 

Table 2 


System geometry, all properties the gas and the solid and operating conditions 


Fluidized bed dimension 


Height 30.00 cm 
Width 10.85 cm 
Number of cell x cell size in x direction 31 x 0.35 cm 
Number of cell x cell size in y direction 20 x 1.5cm 
Fluid properties 
Initially the bed was filled with the N2 
Steam, H20 was fed into the bed. 
Solid properties 
Initially the bed was filled with the carbon nanoparticles 
Agglomerate particle diameter 204 wm 
Agglomerate particle density 440 kg m~? 
Initial solid volume fraction for agglomerate particles 0.44 
Operating conditions 
Pressure latm 
Temperatures 750°C 
Inlet velocity of steam, H2O 3.20cms7! 
Initial bed height 4.50 cm 
Initial weight fraction of C particles 0.95 
Time interval, At 1075 s 


P 21848.52 
Keg = — 6 = 1.5 x 1076 exp (a) (23) 
(Pi) RTs 
Pe = Pen, (24) 
Hz 1.5 x 1076 exp (21848.52/RTs) 


The rate constant, K;,H,, is given in Yoon et al. [27]. However, the 
activation energy used is much lower. Hence, there is essentially 
no methane production. 


Reaction 3: C+ H20 —> CO+H> 
(AHR 31298 k = 32.432 kcal gmol~') (25) 


45, 000 
K: co, = 930 exp (- RTs ) 
(gmol (gmol C)~! atm”! s~!) (26) 


P, 16326.1 
a= Pa O = exp (17.2931 = ) (27) 
* 
(Pho) S 
Po = ee (28) 
2 exp (17.2931 — 16326.1/Ts) 


The numerical values are for Pittsburgh 8 coal. They are similar 
to those of Yoon et al. [27]. Skinner and Smoot [28] give the 
Arrhenius constants for graphite, charcoal and coke. 

Water shift reaction [27]: 


CO + H20 } CO2 + H2 
(AHR 4l208K = —9.838 kcal gmol”') (29) 


—8421.3 
R4 = 0.775 exp (=) po-5—P/250 
G 


Xco, XH 


x (cox ) (1 — We)pcés 


ws 
(gmol cm~? s7!) (30) 


om) (31) 


Kwe = 0.0265 exp (= 
G 


Wen [31] reviewed the details of the shift reaction. 
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5.3. Gasifier fuel cell 


5.3.1. Description of gasifier fuel cell 

A complete derivation of the basic model equation J;= 
(E(P;) — Ve)/ Rete to obtain the current density is given in Dharia 
[32], in Gidaspow’s report [33,34] and his book [35]. 

The fuel cells produce electricity by the electrochemical oxi- 
dation of hydrogen and carbon monoxide. The half-cell reactions 
are: 


e Hydrogen system: 


Anode H,+CO? & H,O+CO, +2e 
Cathode 1/20, +2e° +CO, & COF 
Overall H,+1/20, @H,O 


e Carbon monoxide system: 


Anode CO + COȘ“ ¢ 2C0, + 2e7 
Cathode CO, +1/20,+2e CO7 
Overall CO+1/20, + CO, 


5.3.2. Model 

For the hydrogen-carbon monoxide gasifier fuel cell the con- 
servation of species equations of hydrogen and carbon monoxide 
are added into the gasification code to describe electrochemical 
reactions. The left-hand side term of the continuity equation for 
gas phase does not equal to 0 due to mass changing with gasi- 
fication and electrochemical reactions. They are summarized as 
follows: 


e Continuity equation for gas phase: 


0 
aches) + V - Ogee Ve) 
= Ri(—Mco, + 2Mco) + R:(—2Mp, + Mcu,) 
+R3(Mco + Mp, — Mmo) + Ra(Mco, — Mco 


I 
Mu, Mco (32) 


+My, + M] 
H2 H20) TFA 


e The conservation of species balance in gas phase: 

component balance for H? : 

ð 

gr (Pete tih) + V - (ogee Yn, Vg) 

I 
= Mp, (—2R2 + R3 + R4) — — Mp, (33) 
n Fæ 
component balance for CO : 


0 
g Pees ¥co) + V - (Pgég¥co Vg) 


I 
= Mco Rı + R3 — Ra) = —— Mco 84 
nFa 


The reversible emf of fuel cell, E(Pp,) and E(Pco) are 
obtained from the Nernst equation, calculated as a function of 


concentration of hydrogen and carbon monoxide as follows: 


1/2 
RT Pp, Po 
E(P = Ei — p (V 35 
(Pm) o+ af Pas (V) (35) 
1/2 
RT PcoPo 
E(P, =E ha (y 36 
(Pco) o+ 5p n A (V) (36) 


Rate consumption of hydrogen and carbon monoxide in a fuel 
cell is Ij/njFa (mol cm~? s71), where J; is the current density 
(A cm™?), ni: 2 (number of electron produced per H? mole); 2 
(number of electron produced per CO mole), F: 96,500 C mol7!, 
a: 2mm, thickness of anode channel. 


5.4. Initial condition and boundary condition 


To solve the equations listed in Table 1, appropriate initial 
and boundary conditions are needed. The initial conditions can 
be specified according to the computed system. In this study, we 
simulated the gasifier fuel cells in a batch mode. At the inlet, the 
gas velocity, the gas compositions, gas temperature and pres- 
sure were prescribed. In the reactor, the carbon nanoparticles 
were contained at minimum fluidizing velocity. Therefore, the 
volume fractions, velocities, compositions, temperatures of gas 
and solids phases and pressure were prescribed. The simulations 
were carried out for a two-dimensional fluidized bed. The sys- 
tem geometry, all properties the gas and the solid and operating 
conditions are summarized in Table 2. 


5.4.1. Initial conditions for gasifier fuel cell 

In the simulations, only the anode is simulated. The con- 
centration of oxygen, Po,, is assumed to be constant. The gases 
produced from the fuel cell, carbon dioxide and steam, Pco, and 
Py,0, are also assumed to be constant with no feedback to 
gasification. The mole fractions of oxygen, carbon dioxide 
and steam are 0.21, 0.01 and 0.01, respectively. The standard 
emf of a fuel cell is Eg = 1.2605 — 0.000257. The terminal cell 
potential, Vin, was varied to be 0.8 and 0.6 V. The effective cell 
resistance for a molten carbonate fuel cell, Reff, estimated from 
the slope of cell potential and current density is 1.987 ohm cm? 
[32]. Chuck [13] had shown that the polarization of the carbon 
monoxide fuel cell at 875 °C is the same as that of the hydrogen 
fuel cell at 800°C. Hence here we used the same resistance for 
the carbon monoxide fuel cell. 


5.4.2. Boundary conditions 

At an impenetrable solid wall the gas velocities in two direc- 
tions are generally set to be zero. The no-slip condition was 
applied for solids phase. At the outlet, the mass flux is assumed 
to be continuous. 


6. CFD simulations 


During the energy crisis in 1970s, the design of gasifiers and 
fluidized bed combustors could not be done from first principles 
due to a lack of understanding of gas-particle flow. With funding 
from National Science Foundation and Department of Energy 


D. Gidaspow, V. Jiradilok / Journal of Power Sources 166 (2007) 400-410 


in 1980s and 1990s and the development of fast computers it 
became possible to solve the coupled Navier-Stokes equations 
using the computational fluid dynamics (CFD) approach. 

An example of the acceptance of CFD in the fluidization 
community is the plenary lecture by Kuippers et al. [36] at Flu- 
idization IX. Two plenary lectures in the Flour—Daniel AIChE 
lecture series have summarized the progress made in the last 
decade. Jackson’s [37] 1994 lecture discussed the equations of 
motion and hydrodynamic instability, while Arastoopoutr’s [38] 
1999 lecture emphasized their numerical solution. Commer- 
cial codes, such as Fluent, allow the solution of transient multi 
dimensional problems. Gidaspow [35] had reviewed the theory. 

Table 1 summarizes the conservations of mass, momentum 
and energy equations for the gases and the particles. The con- 
servations of species equations are shown in the gasification 
and gasifier fuel cell sections. The gasification reactions are 
similar to those for a moving-bed coal gasification reactor of 
Yoon et al. [27]. The molten carbonate fuel cell model is based 
on the research at The Institute of Gas Technology and Fuel 
Cell Energy Corporation as reviewed by Gidaspow [34]. In the 
present simulations the gasification kinetics were limited by the 
rates of chemical reactions determined for large coal particles. 
For nanoparticles we expect the rates of reactions to be much 
higher due to their large surface area. 

Figs. 3—7 show our simulations for a fuel cell reactor at 
isothermal conditions, 30.0 cm x 10.9 cm, initially filled with 
close to pure carbon up to the height of 4.5 cm fluidized with 
nitrogen. Fig. 3a shows the carbon distribution after 20 s of reac- 
tions. Fig. 3b and c shows the hydrogen concentration and the 
current density, respectively. In this simulation the inlet steam 
velocity was low, 3.2 cm s~! and the bed was only partially filled 
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Fig. 4. Gasification axial profiles with no current production, 1023 K, 3.2 m s7! 
steam. 


with carbon. Hence, the reaction and current generation were all 
at the bottom of the reactor. The current density in Fig. 3c is 
with no carbon monoxide reaction. In a later section, it is shown 
that the current density is higher when carbon monoxide reacted 
electrochemically. 

Fig. 4 shows gasification with and without the shift reac- 
tion. With the water shift reaction the hydrogen concentration 
increases from 40 to 50%. Similar to the Lurgi gasifier [27] there 
is appreciable carbon monoxide formation. 

Fig. 5 shows the compositions in the gasifier fuel cell fed 
with a low steam velocity, 3.2 cm s7! and a high steam velocity, 
32cms7!. At the low velocity no hydrogen leaves the reactor. 
The product is almost all carbon monoxide. At the high veloc- 
ity, about 40% hydrogen and 50% carbon monoxide leave the 
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Fig. 3. Computed (a) solid volume fraction of carbon, (b) weight fraction of hydrogen and (c) current density at 0.8 V with no electrochemical reaction of carbon 


monoxide, 1023 K. 
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Fig. 5. Concentration axial profiles of a gasifier fuel cell, with no electrochemical 
reaction of carbon monoxide, at high and low steam velocities, 1023 K, 0.8 V. 
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Fig. 6. Current density axial profiles of gasifier hydrogen and carbon monoxide 
fuel cells, operating at high and low steam velocities, 1023 K, 0.8 V. 


cell. In the first study, we assumed that carbon monoxide does 
not contribute to the formation of current. In the molten carbon- 
ate fuel cell it was shown by Chuck [13] that carbon monoxide 
can be oxidized electrochemically to produce reasonable cur- 
rent densities. Morita et al. [12] computed the electrochemical 
work output from the molten carbonate fuel cell by considering 
both the electrochemical oxidization of hydrogen and carbon 
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Fig. 7. Carbon conversion in a batch mode as a function of operating time, at 
1023 K, 3.2cms7! steam. 
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Fig. 8. Heat generation in the fuel cell with hydrogen oxidation only and gasi- 
fication at 0.8 and 0.6 V, 3.2ms7!. 


monoxide. Hence, in the next modeling study we included the 
production of current due to the electrochemical oxidation of 
carbon monoxide. 

Fig. 6 shows the current densities for the low and high steam 
velocities with an electrochemical production of current by 
hydrogen and carbon monxide. At the high steam velocity the 
current density is uniform, except at the inlet. The current den- 
sity with electrochemical carbon monoxide reaction is double 
that with hydrogen oxidation only. 

Fig. 7 shows the weight fraction of carbon conversion as a 
function of operating time in the gasifier only and the gasifier 
fuel cell operating at 0.8 and 0.6 V. Gasification speeds up with 
a fuel cell due to the shift in equilibrium. We plan to simulate the 
gasifier fuel cells in a batch and continuous carbon feed modes 
and find optimum operating conditions. 

Figs. 8-11 show the operation of gasifier fuel cells at 
adiabatic conditions. The initial temperature was 1023 K. 
Figs. 8 and 9 show that we have to supply heat at the inlet of the 
gasifier fuel cell and remove it in the second part of the reactor. 
The corresponding temperature profiles in Fig. 10 shows a drop 
in the temperature near the inlet followed by a rise. The drop in 
temperature in the upper portion of the system is due to the tran- 
sient operation at the low steam velocity. There is no reaction in 
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Fig. 9. Heat generation in the hydrogen—carbon monoxide gasifier fuel cell at 
0.8 V, 3.2m s7!. 
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Fig. 10. Temperature and total heat profiles of hydrogen—carbon monoxide 
gasifier fuel cell at 0.8 V, 3.2ms7!, 
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Fig. 11. Concentration axial profiles of a hydrogen—carbon monoxide gasifier 
fuel cell, at 0.8 V, 3.2ms7! at adiabatic conditions. 


the top part. An optimum steam velocity needs to be determined. 
We also plan to add the cathode stream and the boiler to produce 
steam. 

Fig. 11 shows the concentration profiles at adiabatic condi- 
tions for the gasifier carbon monoxide—hydrogen fuel cell. The 
computed temperature distributions for adiabatic operation show 
an initial drop in temperature due to gasification, followed by a 
rise which will have to be balanced by staggering the cells or by 
other means using the CFD design method. Unlike at isothermal 
conditions shown in Fig. 5, here the carbon monoxide decreased 
to zero at the outlet. Hence, here we have a complete conversion 
of hydrogen and carbon monoxide. 


7. Conclusions 


(1) Traditionally [31] coal gasification processes were opti- 
mized using measured rates of gasification reactions and 
macroscopic type balances of species and energy, such as 
now found in commercial codes like ASPEN. The devel- 
opment of multiphase flow theory [35] and fast computers 
now permits the detailed design of such reactors using com- 
putational fluid dynamics. This study provides an example 


of such an approach using a nanoparticle size CFD model 
of Jiradilok et al. [17]. 

(2) Integration of gasification and molten carbonate fuel cells 
has been previously suggested in the literature. Here, we 
show using CFD calculations that gasification and electric- 
ity production can be done in the same unit, similarly to 
the production of electricity from natural gas with internal 
reforming using the molten carbonate fuel cell. 

(3) Fig. 11 shows that carbon monoxide will react completely in 
the fuel cell with electrochemical reaction. With no reaction, 
Fig. 5 shows that it will leave the fuel cell gasifier partially 
unreacted. Hence, to take full advantage of the new concept, 
carbon monoxide will have to react electrochemically in the 
fuel cell gasifier. 
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